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Introduction 

The  original  research  under  this  contract  was  directed  toward  comparison  of  ocean  bot¬ 
tom  seismometer  and  hydrophone  data  with  ocean  si  b-bottom  seismometer  data  and  with 
appropriate  synthetics  to  determine  the  effect  of  sub-boitom  structure  on  the  signals  received; 
expected  propagation  loss;  predominant  modes  of  propagation;  instrument  signal  distortion  and 
excess  noise;  expected  spacial  coherence;  and  expected  improvement  from  polarization  and 
array  processing.  During  the  first  and  subsequent  years  of  the  contract,  research  emphasis  was 
re-directed  toward  analysis  and  interpretation  of  data  from  the  Columbia-Point  Arena  ocean 
bottom  seismic  station  (OBSS)  emphasizing  the  objectives  of  the  ONR  Accelerated  Research 
Initiative  in  ULF/VLF  Ocean  Seismo- Acoustics.  (  )  «aJ- - 

The  most  extensive  existing  set  of  ocean  bottom  data  on  ULF/VLF  seismic  background 
noise  and  signals  was  obtained  from  OBSS,  which  operated  for  over  six  years  (1966-1972) 
about  200  kin  west  of  San  Francisco  at  a  depth  of  3900  m.  The  OBSS  included  a  Lamont 
long-period  (LP)  triaxial  seismometer  (15  sec  natural  period,  originally  developed  for  lunar 
use);  three-component  short-period  (SP)  system  (1  sec  natural  period);  lo  lg-period  (crystal) 
hydrophone;  short-period  (coil-magnet)  hydrophone;  ultralong-period  (Vibratron)  pressure 
transducer,  thermometer,  current  amplitude  sensor,  and  a  current  direction  sensor.  It  was  con¬ 
nected  to  shore  via  cable  and  LP  and  SP  components  were  recorded  separately  on  FM  tape  as 
well  as  on  seismograph  drum  recorders.  During  this  research,  original  -M  tapes  from  the 
OBSS  long-  and  short-period  (LP/SP)  hydrophones  and  LP  and  SP  three-component  seismom¬ 
eters  were  digitized  to  obtain  ULF/VLF  spectra  and  covariances  during  quiet  and  noisy  times 
and  during  passage  of  vessels  and  earthquake  wavetrains.  These  were  compared  with  available 
source  information  on  sea/swell  condition,  tidal  currents,  earthquakes,  and  shipping.  Analyses 
were  conducted  in  the  frequency  band  .002  to  5  Hz. 

Coherent  energy  peaks  near  0.14,  an  !  0.06  Hz  suggest  fundamental  mode  Rayleigh  wave 
motion  propagating  shoreward.  Coherent  energy  near  .30  Hz  appears  to  be  variable.  Pressure 
variations  near  0.01  Hz  and  lower  frequency  co  Tel  ate  with  wave  heights  along  the  California 
coast  and  appear  to  produce  forced  deformation  of  the  bottom.  Two  intense  NE  Pacific 
storms  with  hurricane  force  winds  occurred  during  one  two  week  time  period.  Time  varia¬ 
tions  of  spectra  and  of  amplitude  and  phase  coherencies  of  the  four-component  OBSS  data 
were  related  to  the  storm  histories  and  to  ocal  weather/wave  conditions  and  used  to  identify 
motion  (seismic  wave)  types  and  directions  of  propagation. 

The  research  supported  by  this  contract  is  summarized  in  the  reports  and  publications 
listed  below.  A  reprint  and  a  preprint  of  two  publications  are  included  as  appendices  A  - 
'  "Particle  Motion  and  Pressure  Relationships  of  Ocean  Bottom  Noise:  3900  m  Depth;  0.003  to 
5  Hz"  and  B  -  "Ocean-Bottom  ULF  Seismo-Acoustic  Noise:  .002  to  .4  Hz." 
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PARTICLE  MOTION  AND  PRESSURE  RELATIONSHIPS  OF  OCEAN  BOTTOM  NOISE: 

3900  M  DEPTH;  0.003  to  5  HZ 

Noel  Barstow,  George  H.  Sutton,  and  Jerry  A.  Carter 

Rondout  Associates,  Incorporated 


Abstract.  Samples  of  ocean  bottom  noise  in  the  frequency 
band  0.003  to  5  Hz  are  analyzed  for  coherency  and  amplitude 
and  phase  relationships  among  pressure  and  the  three  com¬ 
ponents  of  particle  motion.  Data  available  from  the 
Columbia-Point  Arena  ocean  bottom  seismic  station  (38 ' 
09.2 'N,  124'  54.4 'W)  provide  examples  of  different  noise 
conditions.  Coherent  energy  peaks  near  0.14,  and  0.06  Hz 
suggest  fundamental  mode  Rayleigh  wave  motion  propagat¬ 
ing  shoreward.  Coherent  energy  near  .30  Hz  appears  to  be 
variable.  Pressure  variations  near  0.01  Hz  and  lower  fre¬ 
quency  correlate  with  wave  heights  along  the  California  coast 
and  appear  to  produce  forced  deformation  of  the  bottom. 

Introduction 

Possibly  the  most  extensive  set  of  ocean  bottom  data  on 
long  period  seismic  background  noise  and  signals  was 
obtained  from  the  Columbia-Point  Arena  Ocean  Bottom 
Seismic  Station  (OBSS).  OBSS  operated  for  over  six  years 
and  a  number  of  papers  [e.g.,  Auld  et  al.,  1969:  Latham  and 
Nowroozi,  1968;  Piermattei  and  Nowroozi,  1969;  Sutton  et 
al.,  1965:  Nowroozi  et  al.,  1968  and  1969)  have  been  published 
using  OBSS  data  on,  e.g.,  gravity  and  pressure  tides  and  tidal 
currents;  ocean  bottom  microseisms;  Rayleigh  waves  from 
pure  oceanic  paths  and  across  the  continental  margin;  and 
local  and  regional  earthquakes. 

The  Columbia-Point  Arena  OBSS  was  installed  on  18  May, 
1966  at  38  '  09.2  'N,  124  "  54.4  '  W  about  200  km  west  of  San 
Francisco  at  a  depth  of  3903  meters.  It  was  in  continuous 
operation  for  more  than  six  years,  until  11  September  1972. 
The  OBSS  included  a  Lamont  long-period  (LP)  triaxial 
seismometer  (15  sec  natural  period,  originally  developed  for 
lunar  use);  three-component  short-period  (SP)  system  (1  sec 
natural  period);  long-period  (crystal)  hydrophone;  short- 
period  (coil-magnet)  hydrophone;  ultralong-period  (Vibra- 
tron)  pressure  transducer;  thermometer;  current  amplitude 
sensor:  and  a  current  direction  sensor. 


Data 

The  data,  recorded  on  two  7-channel  instrumentation- 
quality  FM  tape  recorders  and  on  seismograph-type  drum 
recorders  and  strip  chart  recorders,  are  currently  stored  at 
Lamont-Doherty  Geological  Observatory.  Short-  and  long- 
period  response  curves  are  published  in  Sutton  et  al.  (1965). 
The  long-period  data  are  digitized  at  8  sps  and  the  short- 
period  data  at  80  sps;  a  four-pole  Bessel  anti-alias  filter  was 
applied  at  a  frequency  equal  to  one-fifth  the  sample  rate.  For 
each  sample  of  signal  or  noise  digitized,  we  also  digitized 
tU.  ■" — meter  calibration  pulses  for  that  date.  After  system 
responses  are  removed  from  simultaneous  noise  samples  of 
long-  and  short-period  data,  for  each  of  tue  three  components 
of  ground  motion,  the  LP  and  SP  noise  spectra  are  very  well 
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matched  between  about  0.1  and  0.5  Hz  (system  noise  limits 
the  spectra  above  0.5  Hz  for  LP  data  and  below  0.1  Hz  for  SP 
data). 

To  further  check  the  accuracy  of  the  data,  we  digitized  a 
few  earthquakes  to  verify  polarities  for  all  eight  components, 
to  compare  mb  magnitude  determinations  with  those  pub¬ 
lished,  and  to  compare  pressure-vertical  relationships  with 
theoretical  values.  The  only  potential  problem  we  encoun¬ 
tered  was  with  the  pressure  data.  Unlike  the  seismometers, 
the  hydrophone  channels  have  no  daily  calibration  pulses. 
Comparison  of  observed  and  theoretical  P^Z  amplitude  and 
phase  relationships  for  digitized  earthquake  P-waves  and  Ray¬ 
leigh  waves  has  convinced  us  that  the  6  dB /octave  high-pass 
corner  at  3  Hz,  included  in  the  original  calibration  curve  for 
the  coil  (SP)  hydrophone,  does  not  exist  and  should  be 
replaced  with  a  0.3  Hz,  6dB/octave  high-pass  corner.  The 
original  3  Hz  corner  represents  the  effect  of  an  hydraulic  relief 
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Fig.  1.  (A)  Instrument-corrected  noise  power  spectra  of  pres¬ 
sure  and  vertical  velocity.  T  he  pressuie  Lola  may  be  ud  to  6 
dB  too  low.  (B)  Vertical  velocity  power  spectra  for  OBSS 
compared  with  spectra  for  continental  stations  under  quiet 
and  noisy  conditions  ('cn  P«.rercon>  t'crc  unpi.-b’irheA 
data;,  i  ne  oilSS  data  are  simultaneous  samples  from  LP 
(solid  line)  and  SP  (dashed  line)  seismometers. 
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Fig.  2.  (A)  Instrument-corrected  noise  power  spectra  of  pres¬ 
sure  (P),  vertical  (Z)  and  horizontal  (Hi  )  motion.  Units  are 
Pa  for  P  and  n m  for  Z  and  H.  (B)  Coherency  spectra,  Z  and 
P  components,  amplitude  and  phase  of  the  data  in  A.;  phase 


convention:  P  leads  Z.  (C)  Same  as  2B.  but  for  Z  and  H 1 
components;  +HI  is  azimuth  246";  phase  convention:  HI 
leads  Z. 


system  required  for  installation  in  deep  water.  Corrected 
pressure  spectra  from  the  LP  and  SP  (with  the  3  Hz  corner 
shifted  to  0.3  Hz)  hydrophones  agree  for  phase  and  shape  of 
amplitude  between  about  0.1  and  0.5  Hz  (the  overlapping  fre¬ 
quency  band),  but  the  SP  hydrophone  data  is  about  6  dB 
higher  than  the  LP  hydrophone  data;  thus,  the  LP  pressure 
power  spectra  presented  in  this  paper  may  be  up  to  6  dB  too 
low. 

A  final  comment  about  the  data:  though  one  channel  on 
both  the  LP  and  SP  FM  tapes  measured  compensation  for 
wow  and  flutter  during  initial  recording,  we  did  not  have  the 
equipment  to  electrorically  compensate  during  playback. 
Instead,  we  subtracted  the  digitized  compensation  channel 
from  each  of  the  data  channels  and  thereby  improved  data 
quality.  Overall  system  noise  limits  background  data  above 
about  5  Hz  and  below  about  0.003  Hz.  (Except  that  the 
long-period  horizontals  record  tills  down  to  at  least  0.001 
Hz.) 

Results 

Pressure  and  verti-e1  vrl.iCly  spec  tin,  luiituu!  feu  instru¬ 
ment  response,  are  shown  in  Figure  1A.  These  spectra  are  cal¬ 
culated  for  a  one  hour  period  during  which  the  bottom 


current  velocity  did  not  exceed  2  cm/sec.  The  pressure  spec¬ 
tra  agree  with  those  obtained  by  Cox  et  al.  [  1 984]  and  Webb 
and  Cox  [  1986},  showing  a  rapid  decrease  for  frequencies 
below  the  "double  frequency"  microseism  peak  near  .14  Hz  to 
a  strong  minimum  between  about  0.03  and  0.1  Hz.  At  times, 
the  "single  frequency"  microseismic  noise,  at  about  .06  Hz  in 
this  sample,  is  nearly  absent,  giving  a  deeper  minimum.  The 
pressure  and  vertical  velocity  spectra  are  generally  similar 
except  that  the  P/Z  ratio  increases  significantly  below  about 
.02  Hz.  The  P/Z  ratio  for  seismic  waves  is  expected  to 
decrease  with  decreasing  frequency  below  about  .07  Hz 
(Bradner,  19631  and  appears  to  be  correct  for  the  .06  Hz  peak. 
Below  about  .04  Hz  P/Z  is  not  consistent  with  predictions  for 
seismic  waves.  H.  Bradner  and  M.  Reichle  (personal  com¬ 
munication)  have  obtained  comparable  results  for  the  P  and 
Z  components  of  the  OBSS  between  about  0.05  and  0.2  Hz  for 
a  different  time  period. 

In  Figure  IB  the  power  spectrum  from  Figure  1A  with  SP 
vertical  data  added  is  compared  with  seismic  noise  on  con¬ 
tinents.  Curves  for  continents  represent  noisy  and  quiet  condi¬ 
tions  on  hard  rock.  ISimilar  curves  for  continents  and  islands 
are  found  in  Murphy  and  Savino,  1975.1  Between  .06  and  5 
Hz.  i tn.  bottom  ncisc  level  p'-trzeH  here  is  generally 

less  than,  but  close  to.  the  level  of  noisy  continental  data; 
below  .04  Hz  it  is  well  above  continental  noise.  The  greatest 
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Fig.  3.  (A)  Instrument-corrected  noise  power  spectra  of  pres¬ 
sure  (P),  vertical  (Z)  and  horizontal  (Hi  )  motion.  Units  are 
Pa  for  P  and  pm  for  Z  and  H.  (B)  Coherency  spectra,  Z  and 
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P  components,  amplitude  and  j.  base  of  data  in  A.  (C)  Same 
as  B.  but  for  Z  and  Hi  components.  Phase  conventions 
same  as  Figure  2. 


"real"  difference  occurs  near  .01  Hz  where  the  OBSS  noise  is 
over  32  dB  above  the  "noisy"  continental  curve  and  it  is  pos¬ 
sible  that  the  "real"  difference  continues  to  decrease  below 
.004  Hz.  How  much  improvement  would  be  obtained  by 
better  bottom-package  design;  shallow  or  deep  burial  within 
the  bottom-sediment;  or  rigid  coupling  within  the  basalt  of 
layer  two  are  important  questions  to  consider  for  the  design 
of  future  OBS  systems  [e.g.,  Sutton  and  Duennebier,  1987]. 
The  rise  in  spectral  level  at  frequencies  lower  than  the 
minimum  is  believed  to  be  caused  by  ultra-long  ocean  waves 
and  by  meterological  pressure  disturbances  on  ocean-bottom 
and  land  stations,  respectively  (Webb  and  Cox,  1986;  Murphy 
and  Savino,  1975]. 

Power  spectra  and  coherency  among  three  components  - 
pressure  (P),  vertical  (Z),  and  horizontal  motion  perpendicu¬ 
lar  to  shore  ( H J. )-  are  calculated  for  two  noise  samples,  21 
June  and  4  July  1966  (Figures  2  and  3).  The  two  samples 
illustrate  time-variable  features  of  microseismic  noise.  Strong 
coherency  in  Figure  2  near  .06  Hz  ("single  frequency" 
microseisms)  coincides  with  a  peak  in  the  power  spectra  of  all 
three  components.  Between  vertical  displacement  and  pres¬ 
sure,  coherency  is  strong  from  about  .06  to  .14  Hz  and  also  at 
.30  Hz.  Corrected  phase  relations  and  amplitude  ratios  for 
the  spectral  and  coherency  peaks  near  .06,  .14,  and  .30  Hz  are 
appropriate  for  fundaiii—.al  mode  Rayleigh  waves. 
Theoretical  results  for  appropriate  velocity  structures  [Lat¬ 


ham  and  Nowroozi,  1968]  indicate  that  at  .14  Hz  fundamental 
mode  Rayleigh  waves  are  near  the  cross-over  from  retrograde 
(at  the  longer  periods)  to  prograde  particle  motion.  Thus, 
the  observed  phase  relationship  between  Z  and  Hi  indicates 
propagation  toward  shore  at  the  OBSS  location  (about  160 
km  offshore)  for  the  "single  frequency"  microseisms. 
Although  Z  -  H  1  amplitude  coherency  is  not  strong  for  the  .14 
and  .30  Hz  peaks,  a  stable  phase  coherency  at  .14  Hz  indicates 
predominant  jr/2  phase  difference,  opposite  the  sign  of  the  .06 
peak.  If  the  .14  Hz  microseisms  are  prograde,  their  propaga¬ 
tion  is  shoreward.  The  phase  difference  between  Z  and  H  1 
around  .30  Hz,  though  messy,  looks  more  like  jt/2  than  -x/2, 
again  suggesting  shoreward  propagation  of  prograde  funda¬ 
mental  Rayleigh  waves.  In  contrast,  the  phase  difference 
between  Z  and  Hi  in  Figure  3  appears  to  be  -jt/2  which 
would  indicate  seaward  propagation.  The  relatively  poor 
coherency  between  Z  and  Hi  for  both  the  .14  and  .30  Hz 
microseisms  suggests  variable  propagation  directions.  Note 
also  in  Figure  3  that  the  "single  frequency"  microscism  is  not 
well  developed. 

In  Figure  4  a  P-Z  coherency  maximum  near  .01  Hz  coin¬ 
cides  with  a  spectral  maximum  most  clearly  observed  from 
the  hydrophone.  The  ratio  of  pressure  to  vertical  velocity  is 
much  too  high  for  fundamental  mode  Rayleigh  waves.  The 
amplitudes  of  these  disturbances  do  not  correlate  in  time  with 
tidal  currents  but  do  correlate  with  ocean  wave  heights 
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Fig.  4.  Instrument-corrected  noise  power  spectra  and 
coherency  spectra  for  Z  and  P  components.  Units  are  Pa  for 
P  and  pm/sec  for  Z.  Coherency  spectra  of  the  same  data  are 
calculated  using  vertical  displacement;  phase  convention:  P 
leads  Z. 


observed  along  the  nearby  California  coast  (Sutton  et  al., 
1965].  The  pressure  signal  is  produced  either  by  long  (shallow 
water)  waves  non-linearly  generated  near  shore  from  the 
ocean  swell,  or  by  differences  in  pressure  beneath  the  high- 
amplitude  vs.  low-amplitude  swell  as  varying  amplitude 
wave-sets  pass  over  OBSS,  or  a  combination  of  the  two.  In 
the  former  caac  shoaling  water  is  required,  whereas  in  the 
latter  case  the  disturbance  would  be  observed  wherever  water 
depth  is  not  large  compared  to  the  "wavelength"  of  the  wave 
sets.  In  either  case  the  velocities  and  wave  lengths  involved 
are  much  smaller  than  for  seismic  signals  of  the  same  period 
and  the  bottom  moves  as  a  forced  deformation. 

Summary 

Results  discussed  is  this  paper  are:  1)  coherent  spectral 
peaks  are  observed  on  vertical  seismometers  and  hydrophones 
near  .01,  .06,  .14,  and  .3  Hz;  2)  the  peak  at  .06  Hz  also  shows 
strong  coherency  with  horizontal  motion  perpendicular  to  the 
coast;  3)  peaks  at  .06,  .14,  and  .30  Hz  have  comparable 
pressure-vertical  velocity  ratios  and  s/2  phase  difference 
appropriate  for  free-running  boundary  waves;  amplitude  and 


phase  relations,  including  horizontal  components,  suggest  the 
.06  and  .14  Hz  peaks  are  predominantly  fundamental  mode 
Rayleigh  waves  propagating  toward  shore,  but  the  .30  Hi 
peak  is  not  as  clear;  4)  a  peak  at  .01  Hz  has  5-10  times 
greater  pressure- vertical  velocity  ratio  than  shorter  period 
peaks;  the  amplitude  of  the  .01  Hz  pressure  peak  correlates 
with  ocean  wave  height  observed  at  nearby  coastal  stations; 
and  5)  the  noise  spectra  observed  from  the  OBSS  seismo¬ 
graphs,  including  the  deep  (variable)  minimum  between 
about  .03  and  .1  Hz,  are  generally  similar  to  published  pres¬ 
sure  spectra  [e.g.,  Webb  and  Cox,  1986]. 
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Abstract 


Observed  spatial  and  temporal  characteristics  of  ULF  ocean-bottom  seismo-acoustic 
ambient  noise  are  required  in  order  to  construct  realistic  quantitative  predictive  models  of  the 
phenomena  involved  Few  such  data  exist  or  have  been  studied,  especially  for  frequencies 
below  about  0.1  Hz.  We  present  analysis  of  noise  data  in  the  band  .002  to  .4  Hz  from  a 
two-week  period,  11/28  to  12/12/67,  recorded  from  long-period,  3-component  seismometers 
and  a  hydrophone  of  the  Columbia-Point  Arena  ocean-bottom  seismic  station  (OBSS,  38° 
09.2'N  -  124°  54.4^,  3903  m  depth).  Two  intense  NE  Pacific  storms  with  hurricane  force 
winds  occurred  during  the  emphasized  time  period.  Time  variations  of  spectra  and  of  ampli¬ 
tude  and  phase  coherencies  of  the  4-component  OBSS  data  are  related  to  the  storm  histories 
and  to  local  weather/wave  conditions  and  are  used  to  identify  motion  (seismic  wave)  types 
and  directions  of  propagation. 

Introduction 

The  general  objective  of  this  research  is  to  improve  our  quantitative  understanding  of  the 
ambient  acoustic/seismic  noise  field  in  the  ocean  as  a  function  of  location  and  time  in  the 
ULF  (0.001  to  1  Hz)/VLF  (1  to  50  Hz)  frequency  band  by  identifying  and  characterizing  the 
noise  sources  and  determining  the  propagation  characteristics  of  the  ambient  noise  field  as 
functions  of  the  ocean  and  sub-bottom  environmental  parameters.  The  specific  research 
reported  here  characterizes  ULF  ocean  bottom  noise  during  two  northeast  Pacific  storms, 
measured  at  a  single  multi-sensor  station. 

The  most  extensive  set  of  ocean  bottom  data  on  long  period  seismic  background  noise 
and  signals  was  obtained  from  the  Columbia-Point  Arena  Ocean  Bottom  Seismic  Station 
(OBSS,  Sutton  el  al.,  1965).  OBSS  was  installed  on  18  May,  1966  at  38°  09.2'N,  124° 
54.4^  about  200  km  west  of  San  Francisco  at  a  depth  of  3903  meters  and  it  was  in  continu¬ 
ous  operation  for  more  than  six  years,  until  11  September  1972.  Though  a  number  of  papers 


(e.g.  Auld  et  al. ,  1969;  Latham  and  Nowroozi,  1968;  Piermattei  and  Nowroozi,  1969;  Sutton 
et  al.,  1965;  Nowroozi  et  al.,  1968  and  1969)  have  been  published  using  OBSS  data,  little 
had  been  published  on  the  background  noise,  especially  outside  the  principal  microseism  band, 
.1  to  .2  Hz,  until  recent  interest  in  ULF/VLF  noise  led  us  to  examine  OBSS  data  specifically 
for  this  information.  OBSS  instruments  primarily  used  for  this  study  are  the  Lamont  long- 
period  (LP)  tnaxial  seismometer  (15  sec  natural  period,  originally  developed  for  lunar  use) 
and  the  long-period  (crystal)  hydrophone.  Records  from  the  three-component  short-period 
(SP)  system  (1  sec  natural  period)  and  short-period  (coil-magnet)  hydrophone  also  were  con¬ 
sulted.  We  digitized  FM  tapes  of  the  long-period  data  at  8  sps. 

We  present  analyses  of  OBSS  seismo-acoustic  ULF  data  for  a  two  week  period,  1 1/28  to 
12/11/67,  during  which  two  NE  Pacific  storms  with  hurricane  force  winds  occurred.  Samples 
of  ocean  bottom  noise  in  the  frequency  band  0.002  to  0.4  Hz  are  analyzed  for  coherency  and 
amplitude  and  phase  relationships  among  the  three  components  of  particle  motion  and  pres¬ 
sure. 


Results 

Figure  1  shows  the  location  and  the  orientation  of  the  horizontal  components  of  OBSS 
and  the  tracks  of  the  two  storms  studied.  Maximum  intensities  of  storm  activity  occurred 
between  about  1200  UT  12/1  and  1800  UT  12/2  for  the  earlier  storm  and  between  about  1200 
UT  12/5  and  0000  UT  12/6  for  the  later  storm  (Mariners’  Weather  Log,  1968).  From  11/30 
through  12/8  (except  12/2)  swell  height  near  OBSS  was  greater  than  15  feet  and  on  11/30  and 
12/3-12/5  combined  sea  and  swell  was  greater  than  20  feet  (historical  conditions  obtained 
from  Pacific  Weather  Analysis).  For  the  period  11/28  to  12/11/67  we  digitized,  somewhat 
arbitrarily,  6-,  12-,  or  4-hour  samples  of  noise  (Figure  2).  These  samples  of  ocean  bottom 
noise  are  used  to  determine  spectra,  coherency,  and  amplitude  and  phase  relationships  among 
the  three  components  of  particle  motion  and  pressure  in  the  frequency  band  .002  to  0.4  Hz. 


Below  about  .003  Hz  some  low  amplitude  samples  may  be  affected  by  system  noise. 

The  pressure  spectra  generally  agree  with  those  obtained  by  Cox  et  al.  (1984)  and  Webb 
and  Cox  (1986)  showing  a  strong  minimum  between  about  0.03  and  0.1  Hz  referred  to  as  the 
"noise  notch";  motion  spectra  from  the  OBSS  seismometers,  especially  the  vertical  sensor,  are 
qualitatively  similar  to  the  pressure  spectra  (Figure  3).  Figure  3  shows  velocity  and  pressure 
spectra  from  three  different  times:  time  1  (0200  UT  11/28)  is  before  the  storms;  time  2  (1200 
UT  12/4)  occurs  after  the  peak  in  double-frequency  (twice  ocean  swell  frequency)  microse¬ 
isms  associated  with  storm  1;  and  time  3  (0200  UT  12/8)  is  near  the  peak  double-frequency 
microseisms  associated  with  storm  2.  It  can  be  seen  that  the  whole  spectral  level  between 
.003  and  .4  Hz  is  raised  for  both  vertical  velocity,  Z,  and  pressure,  P.  For  frequencies  below 
about  .08  Hz  the  horizontal  componer's  are  affected  more  by  local  tidal  currents  (most  likely 
from  direct  turbulent  pressure  on  the  OBSS  package,  Sutton,  1986;  Sutton  and  Duennebier, 
1987  and  1989)  than  by  the  storm  activity. 

The  data  divide  naturally  into  five  frequency  bands  (.004-.02  Hz,  .02-.04  Hz,  .04-.08  Hz, 
.08-. 2  Hz,  and  .2-.4  Hz)  based  upon  observed  characteristics  and  upon  assumed  mode  of  on- 
gin  and/or  propagation.  The  band  divisions  are  labeled  in  Figure  4,  where  amplitude  and 
phase  coherency  between  pressure  and  vertical  motion  are  shown.  Below  we  discuss  each 
band  separately. 

1.  Band  .004-.02  Hz:  In  this  band  it  is  believed  that  the  pressure  variations  result  from 
long  gravity  water  waves  having  wavelengths  greater  than  water  depth  but  much  less  than 
seismic  waves  of  the  same  frequency  (Webb  and  Cox,  1986;  Webb,  1988,  Barstow  et  al., 
1989).  Webb  and  Cox  (1986)  predict  that  this  band  should  be  fairly  constant  with  time;  how¬ 
ever,  we  observe  an  increase  in  P  and  Z  near  .01  Hz  of  about  20  dB  during  the  storms  (Figure 
3).  Sutton  et  al.  (1965)  also  found  a  temporal  correlation  with  ocean  wave  heights  along  the 
nearby  (to  OBSS)  California  coast.  The  pressure  disturbance  should  produce  a  forced  defor¬ 
mation  of  the  bottom  and  near  .01  Hz  we  observe  a  strong  correlation  between  P  and  Z  in 


Figures  4  and  5.  Yamamoto  and  Torii  (1986)  use  such  correlated  pressure  and  motion  in 
shallow  watei  to  estimate  rigidity  vs.  depth  in  the  sub-bottom  as  a  function  of  gravity  wave¬ 
length  (period).  This  theory  assumes  that  variations  are  rlow  enough  that  static  leading  theory 
is  adequate.  The  approximately  3rt/2  phase  between  vertical  displacement  and  pressure  (max¬ 
imum  pressure  in  phase  with  maximum  downward  velocity)  observed  in  Figure  4  indicates 
that  dissipation  (radiation)  is  involved  and  that  the  static  theory  is  not  adequate.  Static  theory 
would  predict  k  phase  (maximum  pressure  in  phase  with  maximum  downward  displacement) 
as  is  also  the  c  tse  for  Rayleigh  waves.  The  pressure-vertical  velocity  .ado  ( P/Z )  near  .01  Hz 
is  24  dB  ±  2.7  S.D.  for  21  samples,  much  too  large  for  fundamenta1  mode  Rayleigh  waves. 
All  P/Z  ratios  were  obtained  by  graphical  measurement  of  individual  power  spectra  calculated 
from  one-huiir  samples  divided  into  8  Hanning  windows  with  62.5  %  overlap. 

2.  Band  .02- ,04  Hz:  In  this  band  we  observe  poor  P-Z  coherence  (Figures  4  and  5). 

The  P/Z  ratio  near  .02  Hz  is  13.2  dB  ±  6.2  S.D.  for  21  samples.  This  avera^**  is  well  above 

that  expected  for  seismic  waves.  Near  .04  Hz,  however,  some  samples  have  ratios  appropriate 
for  fundamental  mode  Rayleigh  waves  (LR).  During  the  storm,  amplitudes  near  .04  Hz 
increased  10-20  dB  (Figure  3).  In  this  band,  which  is  partially  within  the  .03-.  1  Hz  "noise 
notch",  some  low  amplitude  samples  may  be  system  noise  limited. 

3.  Band  .04- .08  Hz:  P-Z  coherence  is  variable  in  this  band,  which  is  also  within  the  .03 

-.1  Hz  "noise  notch"  (Figures  4  and  5).  During  the  storm,  amplitudes  near  .06  Hz  increased 

10-15  dB.  The  P/Z  ratio  near  .06  Hz  for  21  samples  is  6.8  dB  ±  3.1  S.D.;  again,  some  low- 
amplitude  samples  may  be  system  noise  limited.  The  highest  amplitude  samples  (Figure  7) 
have  amplitude  and  phase  relationships  consistent  with  fundamental  mode  Rayleigh  waves. 

This  is  the  frequency  band  of  single-frequency  microseisms  and  of  prominent  oceanic 
Rayleigh  waves  from  earthquakes.  The  single-  frequency  microseisms  are  generally  believed 
to  be  produced  from  direct  pressure  of  shoaling  swell  in  coastal  areas  and  have  been  observed 
mosdy  on  land  instruments  (Oliver,  1962;  Oliver  and  Page,  1963;  Haubrich  and  McCamy, 


1963;  Cessaro  and  Chan,  1989).  Consistent  with  this  are  Webb  and  Constable’s  (1986) 
interpretation  of  ambiguous  results  from  a  two  element  ocean-bottom  array  as  indicating  Ray¬ 
leigh  propagation  from  near  shore.  Spectra  (peaking  near  .06  Hz)  and  coherencies  observed 
from  June  1966  OBSS  data  by  Barstow  et  al.  (1989),  however,  clearly  indicate  fundamental 
mode  Rayleigh  waves  propagating  shoreward  rather  than  seaward  as  would  be  expected  from 
generation  near  the  California  coast. 

Webb  and  Cox  (1986)  noted  the  transient  nature  of  spectral  peaks  for  ocean-bottom 
single-frequency  microseisms  compared  to  more  stable  double-frequency  microseisms.  We 
observe  generally  strong  but  variable  P-Z  coherency  during  the  11/28  -  12/11/67  time  period, 
but  there  are  few  clear  spectral  peaks  in  this  bandwidth.  However,  within  the  three-hour 
period  of  high  P-Z  coherency  near  .06  Hz  (0200-0500  UT  12/8,  Figure  5),  we  are  able  to 
observe  single-frequency  microseisms  by  narrow  band-pass  filtering  the  seismograms  (Figure 
6).  These  microseisms,  peaked  at  .05-.06  Hz,  have  proper  amplitude  and  phase  relationships 
for  LR  and  indicate  shoreward  (from  260°)  propagation,  as  did  our  earlier  results  (Barstow  et 
al. ,  1989).  Also,  we  do  not  observe  a  precise  2  to  1  frequency  relationship  with  the  observed 
double-frequency  microseisms  as  would  be  expected  if  the  same  swell  produced  both  types. 
Oliver  and  Page  (1963)  and  Haubrich  and  McCamy  (1963),  for  example,  at  land  stations 
observed  a  quite  accurate  2  to  1  relationship. 

4.  Band  .08-.2  Hz:  This  is  the  band  of  the  prominent  double-frequency  microseisms. 
We  observe  high  P-Z  coherency  and  a  large  increase  in  amplitude  (30-40  dB)  and  an  increase 
in  period  during  the  storms  (Figures  2  and  3).  The  P/Z  ratio  near  .1  Hz  is  15.2  dB  ±  .8  S.D. 
for  21  samples.  The  amplitude  and  phase  relationships  in  this  band  indicate  fundamental 
mode  Rayleigh  waves  propagating  shoreward  assuming  prograde  motion  based  on  an 
appropriate  velocity  structure  (Latham  and  Nowroozi,  1968;  Bradner  and  Latham,  1972).  The 
evidence,  least  ambiguous  for  double-frequency  microseisms,  is  as  follows.  (The  same  kinds 
of  data  and  reasoning  were  used  to  characterize  single-frequency  microseisms  and  will  be 


discussed  for  wind-wave  microseisms.)  The  7t  phase  difference  between  pressure  and  vertical 
displacement  is  correct  for  Rayleigh  wave  motion  and,  though  shown  only  in  Figure  4,  we 
have  observed  it  in  ail  the  coherence  spectra  between  roughly  .04  and  .4  Hz  (except  at  .08 
and  .2  Hz).  The  P/Z  amplitude  ratios  from  .04-.  15  Hz  are  consistent  with  fundamental  mode 
Rayleigh  waves  for  several  reasonable  velocity  models  (Figure  7).  The  theoretical  points  on 
the  P/Z  ratio  plots  in  Figure  7  are  from  models  5  and  8  of  Bradner,  1963.  These  models  are 
probably  the  closest  of  those  shown  in  Figure  7  to  the  structure  at  OBSS  (e.g.  Sutton  et  al. , 
1971);  model  values  are  given  in  the  table  below. 


Bradner  (1963)  Models  5  and  8 


Model# 

HFm 

p,g/cm3 

Vs, km/sec 

Vp,km/sec 

5 

6 

1.03 

1.52 

1 

1.03 

0.5 

1.52 

6 

2.75 

3.98 

6.9 

CO 

3.09 

4.68 

8.1 

8 

6 

1.03 

1.52 

0.5 

1.03 

0.25 

1.52 

1.5 

2.54 

2.77 

4.8 

5 

2.75 

3.98 

6.9 

CO 

3.09 

4.68 

8.1 

The  H/Z  amplitude  ratios  are  also  consistent  with  fundamental  mode  Rayleigh  waves 
(Figure  8).  Velocity  models  appropriate  for  the  OBSS  site  indicate  prograde  particle  motion 
for  LR  at  frequencies  higher  than  about  .1  Hz.  Assuming  prograde  motion,  then,  the  double¬ 
frequency  microseisms  are  propagating  shoreward  (Figure  9).  In  Figure  9  the  vertical  to  hor¬ 
izontal  coherencies  and  transfer  functions  are  used  to  determine  the  propagation  direction  of 
Rayleigh  waves.  The  phase  coherencies  determine  the  appropriate  quadrant  and  the  ratio  of 
ihe  transfer  functions  determines  the  angle  relative  from  Hll  (Figure  1).  Between  .09  and  .15 
Hz  (double  frequency)  propagation  is  from  the  west  (275°).  In  the  same  manner,  we  deter¬ 
mined  a  propagation  direction  of  260°  for  single-frequency  microseisms  near  .06  Hz  where 
fundamental  mode  LR  is  retrograde. 


The  relationship  between  the  storms  and  the  timing  of  changes  in  double-frequency 
microseisms  is  a  clue  to  microseism  origin.  Using  the  distance  to  OBSS  of  the  closest 
approach  of  hurricane  force  winds  generated  in  storm  1  (Figure  1),  the  calculated  arrival  time 
for  17  sec  water  waves  is  1200  UT  on  12/3.  Similarly  for  storm  2,  the  calculated  arrival  time 
is  2000  UT  on  12/7.  The  earlier  of  two  observed  maxima  for  8  1/2  sec  microseisms  occurs  at 
about  1200  UT  12/3  (Figure  2)  as  predicted  for  generation  near  OBSS  from  interfering  17  sec 
swell  of  storm  1.  The  later-occurring  maximum  period  and  amplitude  of  double-  frequency 
microseisms  is  observed  early  on  12/8,  somewhat  later  than  predicted  for  swell  to  reach  OBSS 
from  the  distances  of  hurricane  force  winds  associated  with  storm  2.  It  is  important  to  our 
understanding  of  the  generation  of  double-frequency  microseisms  that  the  observed  high- 
amplitude,  longest-period  arrivals  from  both  storms  are  later  than  expected  for  Rayleigh  waves 
generated  either  near  the  storm  center  or  at  the  coastal  areas  nearest  the  storm  center.  Thus 
the  time  delay  relative  to  the  storm  activity  suggests  double-frequency  microseism  generation 
near  OBSS.  We  do  not  have  a  simple  explanation  for  the  two-to-three-day  cycle  of  double¬ 
frequency  amplitude  variations  in  Figure  2.  From  11/30  to  12/8  the  low-amplitude  values 
steadily  increase,  and  then  on  12/10  drop  back  to  pre-storm  low  levels.  The  pattern  of  high- 
amplitude  values  is  puzzling,  yet  the  two  longest-period  peaks  in  amplitude  sensibly  relate  to 
the  storm  tracks. 

There  have  been  many  investigations,  mostly  involving  only  land  observations,  of 
double-frequency  microseisms.  They  are  generally  believed  to  be  generated  by  non-linear 
wave-wave  interaction  either  near  a  storm  center  or  upon  reflection  from  a  coastline.  Our 
ocean  bottom  data  seem  to  require  interaction  offshore  and  away  from  the  storm  center. 
Perhaps  the  double-frequency  microseisms  from  both  storms  are  a  result  of  interference 
between  direct  arriving  storm  swell  and  reflections  off  northern  California  and  Oregon  (Figure 
1). 

5.  Band  .2-.4  Hz:  Figure  5  shows  a  stable  P-Z  coherency  maximum  near  .3  Hz  and  a 


minimum  near  .2  Hz  separating  it  from  the  main  double-frequency  maximum.  Amplitudes  at 
.3  Hz  varied  up  to  20  dB  during  the  storms.  Pressure  and  velocity  amplitude  fluctuations  in 
this  band  correlate  well  with  local  wind  and  weather  conditions;  the  generation  of  these 
microseisms  is  related  to  wind  waves.  On  the  original  short-period  photographic  records,  we 
observed  two  pronounced  amplitude  maxima,  12/3  and  12/7.  These  coincide  with  the  two 
maxima  in  local  winds  identified  by  Pacific  Weather  Analysis  for  the  two  week  period  stu- 
died:  30  knots  and  above  from  the  southeast  and  the  northwest,  respectively.  The  P/Z  ratio 
near  .3  Hz  is  6.9  dB  ±  1.0  S.D.  for  16  samples.  We  note  that  observed  P/Z  ratios  above  about 
.15  Hz  do  not  match  fundamental  mode  theory  well  (Figure  7).  Higher  modes  change  ampli¬ 
tude  ratio  and  particle  motion  over  narrow  bandwidths  and  are  highly  sensitive  to  velocity 
structure,  making  it  difficult  to  assess  their  propagation  direction.  If  particle  motion  is  pro¬ 
grade,  propagation  is  seaward,  i.e.  from  Figure  9,  near  .33  Hz  propagation  is  from  the  east 
(100°).  Relatively  poor  coherency  between  vertical  and  horizontal  ground  motion,  however, 
suggests  variable  propagation  direction. 

Summary 

The  pressure  spectra  are  similar  to  those  obtained  by  Cox  et  al.  (1984)  and  Webb  and 
Cox  (1986)  showing  two  maxima  separated  by  a  strong  minimum  between  about  0.03  and  0.1 
Hz;  motion  spectra  from  the  OBSS  seismometers,  especially  the  vertical  sensor,  are  similar  to 
the  pressure  spectra.  Coherent  energy  from  0.1-0.15  Hz  and  intermittently  near  0.06  Hz  indi¬ 
cate  fundamental  mode  Rayleigh  wave  motion  propagating  shoreward.  The  timing  of  high 
amplitude,  longer-period  double-frequency  microseisms  is  consistent  with  microseism  genera¬ 
tion  near  OBSS  rather  than  closer  to  the  storm  center.  Background  noise  levels  around  .06 
Hz,  are  intermittently  coherent  between  P  and  Z  during  the  two  week  period  of  this  study,  but 
rarely  do  they  produce  spectral  peaks.  Coherent  energy  near  .3  Hz  is  consistently  high  and  the 
highest  amplitudes  correspond  to  the  highest  local  winds  for  the  period.  The  P/Z  ratios  show 
that  these  wind-wave  microseisms  are  not  fundamental  mode  Rayleigh  waves. 


Coherent  energy  near  .01  Hz  varies  with  time,  but  is  well  above  random  coherency 
(estimated  from  coherency  between  samples  of  P  and  Z  taken  at  separate  times).  Another 
result  from  both  our  previous  study  of  fairly  quiet  conditions  (Barstow  et  al. ,  1989)  and  the 
current  study  of  severe  storm  conditions  is  a  90°  phase  shift  near  .01  Hz  between  pressure  and 
vertical  displacement,  indicating  dissipation  or  radiation,  in  the  forced  deformation  of  the  bot¬ 
tom  from  long  gravity  water  waves. 
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Location  and  orientation  of  OBSS  and  tracks  of  two  1967  hurricanes.  Thickened 
lines  show  duration  and  location  along  the  storm  tracks  of  hurricane  force  winds. 
Storm  1  (11/29-12/5);  closest  approach  of  center  of  hurricane  winds,  A=15.7  . 
Storm  2  (12/3-12/8);  closest  approach  of  hurricane  winds,  A=19.3  . 


Figure  2.  Measured  from  the  original  long-period  photographic,  continuous  analog  records 
are  predominant  period  (circles)  and  average  peak-peak  trace  amplitude  (trian¬ 
gles)  for  vertical  component  microseismic  background  noise  during  times  of 
amplitude  maxima  and  minima  (usually  lasting  several  hrs).  Horizontal  bars 
represent  time  samples  of  background  noise  digitized  from  original  FM  tapes. 
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Figure  3.  Pressure  and  three-component  velocity  power  spectra  from  three  time  periods:  1 
before  (11/28),  2  during  (12/4),  and  3  near  a  peak  (12/8)  in  microseism  storm 
activity  associated  with  two  NE  Pacific  hurricanes  in  December  1967.  Note  that 
the  whole  spectrum  is  raised  for  pressure  and  vertical  motion  and  only  that 
above  .08  Hz  for  horizontal  motion.  Spectra  were  calculated  from  one-hour 
samples  divided  into  8  Hanning  windows  with  62.5%  overlap. 
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Figure  4.  Amplitude  and  phase  coherency  between  vertical  displacement  and  pressure  dur¬ 
ing  peak  of  microseism  storm.  An  estimate  of  random  amplitude  coherency  is 
given  by  the  generally,  lower- amplitude  spectrum  obtained  between  samples  of  P 
and  Z  taken  at  separate  times;  above  .02  Hz  the  dashed  line  connects  the  largest 
maxima  in  the  spectrum.  7t  phase  (maximum  pressure  in  phase  with  maximum 
down  displacement)  is  correct  for  Rayleigh  wave  motion.  3re/2  phase  near  .01 
Hz  (maximum  pressure  in  phase  with  maximum  down  velocity)  is  not  consistent 
with  simple  static  loading  from  long  gravity  water-waves.  Numbered  frequency 
bands  correspond  to  text.  This  and  succeeding  coherency  spectra  were  calculated 
from  3  Hr.  samples  divided  into  16  Hanning  windows  with  62.5%  overlap. 


Figure  5.  Amplitude  coherency  between  pressure  and  vertical  motion  (.002  to  .4  Hz)  dur¬ 
ing  two  week  period  including  two  hurricanes.  Note  the  variable  coherency  in 
the  ‘single  frequency’  microseism  band  between  about  .05  and  .1  Hz  and  the 
consistently  low  coherency  near  0.2  Hz  and  also  from  .02-.04  Hz. 
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Figure  6.  Pressure  and  vertical  velocity  band-pass  filtered  from  .04-.06  Hz  to  en’  ance 
single-frequency  microseisms;  all  traces  start  at  0300  UT  on  12/8.  Amplitude 
scale  applies  only  to  time-expanded  traces,  for  which  pressure  is  plotted  positive 
downward. 
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Coherency  between  vertical  and  horizontal  motion  (13-1600  UT,  12/7)  and  ratio 
of  the  two  vertical  to  horizontal  transfer  functions.  Assuming  prograde  Rayleigh 
motion,  stable  +  or  -  tc/2  phase  identifies  the  proper  quadrant  of  approach  and 
tan-1(Hl/H II)  gives  the  angle  from  Hll. 
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